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We study the two dimensional XY-model with high precision Monte Carlo techniques and investigate the 
continuum approach of the step-scaling function of its finite volume mass gap. The continuum extrapolated results 
are found consistent with analytic predictions for the finite volume energy spectrum based on the equivalence 
with sine-Gordon theory. To come to this conclusion it was essential to use an also predicted form of logarithmic 
decay of lattice artifacts for the extrapolation. 



1. Introduction 

Interacting quantum field theories in four space 
time dimensions are in all relevant cases impossi- 
ble to solve exactly and results are obtained only 
in suitable approximation schemes. The situation 
is different with two dimensional systems which 
often may be (partially) solved exactly. Such a 
system is for instance the sine-Gordon model for 
which detailed predictions for the finite volume 
energy spectrum exist |ll2l3j . These predictions 
can be extended to the two dimensional 0(2) non- 
linear cr-model which is believed to lie in the same 
universality class as the sine-Gordon model at a 
special value of its coupling. For the derivation 
of these exact results one adopts some not rigor- 
ously provable conjectures at intermediate steps, 
therefore a numerical confirmation is desirable. 

We study the massive phase of the XY-model 
with Monte Carlo techniques and extract the fi- 
nite volume mass gap from a time-slice (zero mo- 
mentum) correlation function. The continuum 
extrapolation is performed according to a predic- 
tion by J. Balog 0]. In the continuum limit some 
points of the step-scaling function of the Liischer- 
Weisz- Wolff (LWW) coupling ^5^ are compared to 
values obtained from solutions of the Destri-de 
Vega (DdV) equation. A more detailed account 
on the subject is given in |^. 



2. Sine-Gordon theory in finite volume 

A complete description of the exact spectrum 
of the sine-Gordon model in finite volume is pro- 
vided by the DdV non-linear integral equations. 
Following [7] we solve the equations iteratively 
and refer to |S] for further details of this proce- 
dure. We are interested in the point where the 
sine-Gordon model and the 0{2) non-linear cr- 
model coincide, that is at PsG = VSn. For a 
chosen volume ML (expressed in units the infi- 
nite volume mass gap M) we solve the equations 
for the ground state and for the first excited state 
in order to calculate the finite volume mass gap 



M{L)^Ei{L)-Eo{L) 



and hence the LWW coupling 



f = 2M{L)L . 



(1) 



(2) 



At the doubled volume a point of the step scaling 
function 



(3) 



'presented by T. Korzec 



is obtained. A number of such results is given in 
the second column of table ^ 
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3. Numerical work 

We simulate the two dimensional XY-model 
with standard action 



si 



(4) 



on L/ax5L/a lattices. We apply periodic bound- 
ary conditions in the spatial direction and free 
boundary conditions in the temporal direction. 
As Monte Carlo algorithm we choose the highly 
efficient single cluster algorithm 8 which does 
practically not suffer from critical slowing down. 
To measure the time slice correlation function 

G(r) = {S{t) ■ S{t + t)), S{t) = i ^ s{x, t) , 

X 

(5) 

we employ Hasenbusch's improved estimator ^ . 
From the correlation function we extract the finite 
volume mass gap 

G(r) 



M{L)a 



In 



(6) 



_G(t + 1 

and obtain a value of the LWW-coupling (f'{L). 
To get one point of the lattice step-scaling func- 
tion 



n'^,f{L),a/L)^f{2L) 



(7) 



we keep (3 fixed and measure the LWW coupling 
on a lattice of doubled size. The procedure is 
repeated several times for the same value of the 
coupling but different lattice resolutions a/L. An 
extrapolation to a/L = yields one point of the 
continuum step-scaling function that can be com- 
pared to the predicted value lO . For that purpose 
the numerical value of ML has to be determined 
first, at which the DdV equations produce the 
same coupling. 

For models with a Kosterlitz-Thoulcss phase 
transition |10lll| like the XY-model, J. Balog has 
predicted the leading lattice artifacts to be uni- 
versal and to decay very slowly, i.e. proportional 
to inverse powers of the logarithm of the infinite 
volume correlation length. Lattice artifacts of the 
step scaling function have the form 



where ^ is the infinite volume correlation length, c 
is universal and can be calculated for each volume 
and L/ is a non-universal constant [U — 1-3(1) 
for the standard action T?). Corrections to this 
formula are of order ln(^)~'*. 

4. Results 

We have performed our calculations at four dif- 
ferent values of the LWW coupling. Table sum- 
marizes our results. The lattice artifacts predic- 
tion includes information about the constant c 
which is also listed in the table. Fig. ^ corre- 
sponds to the second line of the table, plots for 
the other points look similar. The theoretical pre- 
dictions are compatible with continuum extrapo- 
lated lattice results. The small differences be- 
tween cth and CMC may be explained by sublead- 
ing cutoff effects in the Monte Carlo data. The 
knowledge of the form of lattice artifacts was es- 
sential to obtain this result to the precision that 
is reached. 
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Figure 1. Comparison of MC-data with a nu- 
merical solution of the DdV equation at Tp' = 
1.7865. The spatial extents L/a oi the lattices 
were 10, 20, 40, 80, 120 and 160. The smallest 
lattice was discarded in the fit. 
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Table 1 

At different values of the LWW-coupling theoretical predictions (DdV) for the step-scaling function are 
compared with numerical results (MC). Also the slope of the fit c as predicted by theory (th) is compared 
to the numerical value. The last column lists the values of the extrapolations. 



f 


CTDdv(2,5^) 


o-Mc(2,5^) 


Cth 


CMC 


xVdof 


3.0038 


4.3895 


4.40 ±0.02 


2.6176 ±0.0002 


2.4 ±0.6 


2.51/3 


1.7865 


1.8282 


1.829 ±0.007 


5.30 ±0.01 


4.8 ±0.5 


0.73/3 


1.6464 


1.6515 


1.657 ±0.003 


5.4 ±0.2 


4.3 ±0.3 


0.35/3 


1.6020 


1.6029 


1.608 ±0.004 


5.5 ±1.5 


4.4 ±0.5 


0.90/3 



5. Conclusions 

We have investigated the massive scaling limit 
of the XY-model by means of the step-scaling 
function of the LWW coupling and found it to 
be consistent with continuum predictions based 
on the equivalence with sine-Gordon theory. A 
predicted form of lattice artifacts was also con- 
firmed and has been essential to find the agree- 
ment. An extrapolation with powers of a would 
have led to a significantly different continuum re- 
sult in spite of a perfectly reasonable looking fit. 
We would like to recall here that deviations from 
this Symanzik behavior have also been found in 
the asymptotically free 0(3) non-linear cr— model 
where they have not yet been understood jl3114j . 

Our results fit well into the picture that was 
drawn in |12|. where among other things the 
renormalized 4-point coupling and 2-point cor- 
relation functions in the continuum sine-Gordon 
model were compared to their lattice counterparts 
in the XY-model. 
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